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Abstract: We study the rare decays of charged vr mesons, vr"^ — > e~^e~^fj, v^ and vr^ 



^ \ e^^~e^Ve induced by a sterile neutrino A^ with a mass in the range r?!^ < ttitv < it^-k- The 

C^ ' first process violates Lepton Number by two units and so occurs only if N is Majorana, 

ly-s ' while the second process conserves Lepton Number and occurs irrespective of the Majorana 

^^ ' or Dirac character of N. We study a way to distinguish the Majorana vs. Dirac character of 

Cn . N in these processes using the muon spectrum. We also find that the branching ratios could 

^vq . be at the reach of high luminosity experiments like Project X at FNAL or any proposed 
neutrino (or muon) factories worldwide. 
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1 Introduction 

One of the outstanding issues in neutrino physics today is to clarify the Dirac or Majorana 
character of neutrino masses. The discovery of neutrino oscillations indicates that neutrinos 
are massive particles with masses likely to be much smaller than those of charged fermions 
[1-8]. This fact provides an important clue on the existence of a more fundamental physics 
underlying the Standard Model (SM) of particle physics, because neutrinos are naturally 
massless in the SM. 

If neutrinos are Dirac particles, they must have right-handed electroweak singlet com- 
ponents in addition to the known left-handed modes; in such case Lepton Number remains 
as a conserved quantity. Alternatively, if neutrinos are Majorana particles, then a neu- 
trino is indistinguishable from its antiparticle and Lepton Number is violated by two units 
(AL = 2) in some processes that involve neutrinos. The experimental results to date are 
unable to distinguish between these two alternatives. 

There have been several attempts to determine the Majorana nature of neutrinos by 
searching for AL = 2 processes. The most prominent of these processes are neutrino-less 
double beta decays (Oz^/3/3) in nuclei, which are regarded as the most sensitive way to look 
for Lepton Number violation (LNV) [9-12] (for a more recent review see [13]). However, it 
has long been recognized that, even though the Oi^/3/3 experiments are very sensitive, the 
extraction of the neutrino mass scale from nuclear Oz^/3/3 is a difficult task, because reliable 
information on neutrino properties can be inferred only if the nuclear matrix elements for 
Oi//3/3 are also calculated reliably. Even in the most refined treatments, the estimates of 
the nuclear matrix elements remain affected by various large uncertainties [14-19]. 

Another way to detect the Majorana nature of neutrinos is to consider specific scat- 
tering processes with possible lepton number violation in hadronic accelerators [20-23]. 

Yet another avenue to detect the Majorana nature of neutrinos is to study AL = 2 rare 
meson decays [24-29]. In our previous work [29] we have studied rare K, D and B meson 



decays that are sensitive to neutrino masses above m^, up to the masses of B mesons. 
The expected rates for the different decays are quite diverse, due to the kinematics of the 
processes and the current bounds on neutrino masses and mixing parameters. In this paper 
we want to focus our study on the AL = 2 decay of charged vr mesons. This process is 
sensitive to a rather narrow neutrino mass range (between m^ and mvr). However narrow, 
this mass range has good discovery potential in proposed high intensity experiments, such as 
Project X [30] of FNAL or at other neutrino factories or muon factories/cohiders proposed 
worldwide. These kinds of experimental facilities will have very large samples of pion 
decays, much larger than the samples of any other meson decays, so they could be sensitive 
to much smaller decay rates, and so increase the potential for discovery of sterile neutrinos 
in this mass range. 

Specifically, here we consider the decay vr^ — >■ e'^e'^ ^~D^. The observation of this 
exclusive decay will clearly indicate Lepton Number violation. However, because the final 
neutrino flavor cannot be observed in the experiment, a similar but Lepton Number con- 
serving (LNC), AL = (albeit Lepton Flavor violating) process, Tr"^ — >■ e"'~/x~e^fe, cannot 
be easily separated from the signal. We explore how to distinguish these two decays in 
order to confirm the AL = 2 process. We must add that, by specifically choosing /x" in the 
final state of the decay of a vr"^ (i.e. a muon with opposite charge to that of the pion), one 
avoids a serious background coming from the radiative decays vr"'" — > /U"*"!^^ + 7*(— )• e+e"). 

The AL = 2 pion decay we propose can only occur via a Majorana neutrino in the 
intermediate state, and thus its experimental observation can establish the Majorana char- 
acter of neutrinos and the absolute scale of their masses in much the same way as in nuclear 
Oz^/3/3 decays, but there are some essential differences. While the theoretical uncertainties 
in meson decays are much easier to handle than in nuclear Oi^/3/3 decays, the latter are more 
realistic options for experiments, as one can count with macroscopically large samples of 
decaying nuclei. Indeed, AL = 2 meson decay rates, at least in the case of light neutrinos 
(TTijy < 2 eV), are prohibitively small for any beam experiment where the mesons must be 
produced. 

Nevertheless, in the case of a Majorana neutrino (denoted as N) with a mass in the 
intermediate range m^ < rriN < m^, a decay like vr"*" — )• e"'"e"'"/x~i>^ is dominated by a 
resonantly enhanced s-channel amplitude [25, 29, 31]: the intermediate neutrino N goes on 
its mass shell, making the decay rate large enough to be within reach of future experiments. 

Moreover, since the intermediate neutrino N is on its mass shell and N is generally long 
lived, the process actually separates in the sequence vr"*" -^ e'^N followed by A^ — t- e^fi^u. 
The study of the first subprocess may already allow to improve the current upper bounds 
on the magnitude of the lepton mixing element l-BeA^j for mi\f in the specified mass range. 

Besides the effect on decay processes, there is much motivation to search for non- 
standard neutrinos with masses of the order of 0(10*^ — 10^) MeV because their existence 
has nontrivial observable consequences for cosmology and astrophysics, in particular, they 
are presumed to be a component of the Dark Matter in the Universe [32, 33]. Actually, im- 
portant bounds on neutrino masses and their mixing with the standard neutrinos arise from 
Cosmology and Astrophysics [34]. There are also laboratory bounds on sterile neutrinos 
coming from their contribution to various processes that are forbidden in the SM, but those 



bounds turn out to be much weaker than the cosmological and astrophysical bounds, except 
in specific cases where the latter become inapphcable (see [35] and references therein). 

We should remark that neutrinos in our mass range of interest must be sterile with 
respect to weak interactions. From cosmological constraints we know that the standard 
electron neutrino mass is below 2 eV [36], while neutrino oscillation experiments tell us 
that all three neutrino masses differ from one another by much less than that value [37] 
(for a recent review see [38]), and measurements of the Z^ boson width constrain the 
number of standard neutrinos to three. Consequently, all neutrinos with masses above 2 
eV are assumed to be non-standard, and thus sterile with respect to weak interactions. In 
this work we will denote these neutrinos by the letter N, and reserve the letter z/ for the 
standard neutrinos. 

In section 2 we show the calculation of the decay rate of interest, vr"^ -^ e^e^ ii~v^, 
and its related "background" Tr"*" -^ e+e+/u~t'e. We numerically explore the processes and 
discuss the results in section 3, in particular to study the ways to distinguish the Majorana 
vs. Dirac character of the intermediate neutrino using the muon spectrum. In section 4 
we state our conclusions. 

2 Calculation of 7r+ — ^ e~^e~^fi~i'^ and 7r+ — )■ e~^e~^fi~i'e 

The decay 7r+ — )■ e^e^ ii~i>^ is forbidden within the SM as it violates lepton number by 
two units. However, if we assume the existence of a Majorana neutrino A^ with a mass in 
the range m^ < rriM < rriT^ and which mixes with the standard flavors, this decay becomes 
dominated by the amplitude shown in figure 1, where the intermediate Majorana neutrino 
N goes on its mass shell, i.e., it is produced and then decays. Contributions where N is 
not on its mass shell are prohibitively small for any foreseeable experiment, so we will not 
consider them here. Since the final state neutrino flavor is not detected in these kinds of 
experiments, the process vr"*" — t- e'^e^ ^~V(,^ shown in figure 2, should be added in the rate. 
This latter process is also mediated by a neutrino A^, but irrespective of its Majorana or 
Dirac character. Consequently, if A^ is a Majorana neutrino, both processes should occur, 
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Figure 1. Diagram for the Lepton Number violating process vr"* 
Majorana neutrino A'' of mass between m^ and m^. 
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Figure 2. Diagram for the Lepton Flavor violating (yet Lepton Number conserving) process tt"* 
e^e'^ H^Ue mediated by an on-shell neutrino A'' of mass between m^ and rrin- 



while if it is a Dirac neutrino, only the latter process is possible. 

Given that the decay of the intermediate neutrino A^ is necessarily very weak, the two 
positrons in the process can actually be distinguished by their spatially displaced vertices, 
so there is no need to consider the diagram with crossed positron lines (this would not be 
the case if A^ were off-mass shell) . ^ 

We can estimate the expected separation of the vertices by calculating the lifetime 
of the neutrino N. If A^ is a Dirac neutrino, its main decay modes are given in figure 3, 
mediated by (a) charged and (b) neutral weak currents. There are additional (smaller) 
contributions from radiative channels which we neglect. 



(a) (b) 

Figure 3. Main decay of a massive neutrino A*' with tum < ttItt: (a) charged current channels, where the 
charged leptons can be e"'"e~, e'^fi^ or n'^e^ , and (b) possible neutral current channels, where the charged 
lepton pair can only be e'^e~ . 

The dominant charged current mode is A'^ — t- e~^e~u, while the channels with a muon in 
the final state are suppressed by an extra factor /(m,^/?7i^), where f{x) is the well-known 
3-body function that also appears, for example, in muon decay 

f{x) = l-8x + 8x^ -x'^-Ux'^lnx . (2.1) 

This factor, when ttin increases from ?tt,^ to m.,^, decreases from 7.3 x 10"^ to zero, so we 
can neglect the channels with a final muon. Concerning the neutral current channels, these 
appear only in some models of neutrino mixing, and include mixing elements other than 
the Bijy. In models where they appear, they are at most comparable in size to the charged 
ones. We will then consider these neutral channels only implicitly by a factor /C, as shown 
below. 

Now, if A^ is a Majorana neutrino, the channels charge-conjugate to those of figure 
3. a contribute as well, and therefore the witdh r^r is twice as large as in the case of Dirac 
neutrino. Therefore, the decay width of the neutrino A^ in our cases of interest can be 
expressed in general as: 



> ^ r{N -^ e+e-u) CK ^ \BeNV ^ ^ C /C , (2.2) 
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where C = 1 or 2 if A^ is a Dirac or Majorana neutrino, respectively, and /C represents a 
correction to include the neutral current interaction channels shown in figure 3.b (/C would 
be unity if the neutral channels were absent). Using the results of ref. [28], Appendix C, 



^If the crossed amplitude is included for an on-shell A'', the interference term turns out to be zero and 
the two absolute squares give the same contribution, so with the inclusion of the symmetry factor (1/2) 
one obtains the same result as with a single diagram. 



/C turns out to be 

/C = 1.757 + 0.606 ^' ^ ' ' ' ' , 2.3 

\BeNr 

where sin^ 6w = 0.231 was used. The decays contained in this expression are A^ — > i'ie~^e~ , 
N — > i>ie^e~ and N — ?> ViVjUj, where i, j = 1, 2, 3 label the light mass eigenstate neutrinos. 
We denote by B^^ the mixing coefficient between the standard ffavor neutrino v^ {i = 
e, /i, r) and the (heavy) mass eigenstate A'': 

3 

1^1 = Yl ^i^j ^i + ^^^^ ■ (2-4) 

Only one heavy neutrino is included, for simplicity of notation. The extension of our results 
to scenarios with several heavy neutrinos is usually straightforward. 

Numerically, the decay distance CTjsf for a such a Majorana neutrino ranges from about 
300/|i?eArp meters for jnjv = in^^ to about 4 times longer for jnjv = Jtt,^, without considering 
relativistic dilations. Taking into account the current upper bounds |i?eAf P ^ 10"^ [28, 39], 
decay lengths above 10^ meters are expected, so in most cases the vertex separation is 
enormous. For a detector of length L (with L <C ctat), the probability of a neutrino A^ to 
decay inside the detector is Pn = L/jctn (e.g. Pn ~ 10~^/7 for a detector about 100 m 
long), where 7 is the relativistic dilation factor. 

Clearly the search for a massive neutrino N in these processes should first be done by 
looking at the energy spectrum of the positrons produced in the primary vertex, i.e. in a 
search for the pion decay mode n^ — >■ e+A^. For mp^ in the range {m^,mTr), the energy 
of these positrons (in the pion rest frame) is rather low: E^ = {rn^ — mjf)/2m-,r, which is 
at most 30 MeV, for /tt-tv = in^. Notice that, because rriTv is not much less than 771,^, this 
decay is not chirally suppressed (unlike 7r+ — )■ e~^Ue). Using the definition of the pion decay 
constant /vr in terms of the charged axial current, (0|'U7^(1 —J5)d\'7r'^{p)) = i/ttP^, and the 
usual notation for the mixing of the massive neutrino A^ with a lepton flavor i by Bin, the 
decay rate of this mode is: 



1 rf2 f2i,r i2m%{ml 



m 



2 ^2 



r(vr+ ^ e+iV) = -^G^/jf|K, r-^^"'"% '"^^^ \B,n\^ (2.5) 

OVr 777.^ 

which gives a branching fraction (approximating the total width of the pion by the decay 
rate of its dominant channel 7r+ — )> n'^i/^i): 

2 ^2 \2 



mj^[m:;^ — m 
m?{m'^ — m 



Br[.- -. e-N) ^ "^^^^ Z^ l^e^P- (2-6) 
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Consequently, the suppression of this low energy positron peak comes almost entirely from 
the mixing element |-BeArp, except for an additional phase space suppression (1 — m^/rri^)^ 
if TTiTv is near the pion mass threshold. There is no chiral suppression here, unlike in 
TT'^ — )■ e+fe, because A^ is heavy. 

As a preliminary experimental analysis, the study of the positron energy spectrum in 
the primary decay vr^ -^ e^N, specifically the search for bumps in the spectrum below the 



peak coming from the standard process tt'^ — >■ e'^i'e, will allow the discovery of the heavy 
neutrino N, or otherwise to put more stringent upper bounds on the lepton mixing element 

\BeN\- 

Although the observation of the positron in the primary decay vr"^ -^ e^N can lead to 
the discovery of the heavy neutrino, it cannot tell about its Majorana or Dirac character. 
For this latter purpose, we must observe the subsequent N decay. Let us then consider the 
Lepton Number violating decay of interest, tt~^ — )■ e~^e~^fi~v'fj_, depicted in figure 1. In what 
follows we neglect the electron and standard neutrino masses. 

Since the process is dominated by the intermediate neutrino N on mass shell, the 
transition probability can be estimated using the narrow width approximation. The details 
of the derivation of the rates are given in Appendix. After integrating the transition rate 
over all final particles momenta but the muon energy, we obtain the spectrum of the muon 
energy E^ [see eq. (A. 18)]: 
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r(vr+ ^ e+e+fi-i?^) = / dE^ (2.7) 
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where Tj^j is given in Eq. (2.2), with C = 2 as here A^ is a Majorana neutrino. As a reference, 
the general differential decay rate dT/dEi for a general pseudoscalar decay of this type, 
M"*" — >■ £fi2^~ij£, with general nonzero masses of the particles, is given in eq. (A. 7) in the 
Appendix. 

The appearance of r^r in the denominator of these expressions is due to the domi- 
nance of N on its mass shell. Indeed, as in all narrow width approximations, the rate is 
proportional to the branching ratio of the N decay subprocess, namely: 

r(7r+ -^ e+e+fi'u^) = r(7r+ -^ e'+N) ■ Br{N -^ e+fi-Uf,), 

where r(7r^~ — )• e''"A^) is given in eq. (2.5), Br{N — )■ e^ ^~v^) = r(A^ -^ e'^ ^^u^)/Tn, and 
r(A^ -^ e^ li~Ufj) is just like the expression for r(A^ -^ e^e~u) given in eq. (2.2) times a 
factor /{ml/mj^) coming from the phase space with a muon instead of an electron. Given 
the total width r^r in eq. (2.2), the resulting branching ratio is: 

Br{N ^ e+^l-u) = l^f{ml/m%), 

This branching ratio, for m,N in the range {m^^m.-,^) (and C = 2 for a Majorana N) has a 
maximum value -^ 4 x 10~^ for m^j^i = m,.,^, dropping quickly several orders of magnitude 
for lower masses, and vanishing at m,N = m^. On the other hand, the production of A^ in 
the pion decay is phase-space suppressed as m,]\f — > ?tt,^. 

Now, TT~^ — 7> e^e~^fi~v'f^ occurs only if A^ is a Majorana, not Dirac, neutrino. However, 
as previously mentioned, since the final neutrino fiavor is not experimentally detectable, 
there is a "background" for this process that is mediated by a neutrino A^ regardless of its 
Majorana or Dirac character: vr^ -^ e^e^fi'Ve- The diagram for this process is shown in 



figure 2. Following a similar procedure as before, for this process we find the muon energy 
spectrum in the A^ rest frame [see eq. (A. 20)] to be: 



^N+"^fi. 



r(7r+ -^ e+^-e+i/e) = / dE^ (2.8) 

JlUfj, 

32^4 ^3 r^ |("^^ + "^m)^m - ^rnN{2E^ + m^)| ^ El - m2. 

Consequently, when the intermediate neutrino N is Majorana, the measurable muon 
energy spectrum is represented by the sum of expressions (2.7) and (2.8): 



\BeN?E^{m% + m^ - 2mNE^) + \B^n\ 



-Ef,{m% + ml) - -tunEI - -niNml 



where E^ is the muon energy in the neutrino N rest frame, and varies between m^ and 
{rrij^ + m1)/{2m]y). For the decay width Fat, we must use eq. (2.2), with C = 2 because 
here A'^ is Majorana. 

Notice that in this spectrum the differential rate is defined in the frame of the initial 
particle, i.e. the pion, while the variable E^ is the muon energy in the neutrino A^ rest frame. 
If one wanted to express dT also in the A^ rest frame, we should include the relativistic 
time dilation factor I/7 = 2mT^mi\] /{m^ + m?^) for the moving pion in the A^ frame, a 
factor which is between 0.9625 and unity for rriN between m^ and m-,^. To avoid possible 
confusions, we normalize the spectrum relative to the pion decay width, thus obtaining a 
spectral branching ratio: 

— — ^ — (vr^^e^e^/i v) = — ^(vr^ ^ e^e^/x v) 

dEfj_ i ^+ dEfj_ 

,2 \2 
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where we have approximated the total width of the charged pion, F^+, by its dominant 
(by far) decay channel, F(7r''" — ?• /i"*"f^j). Integrating this spectrum over the muon energy 
we obtain the corresponding branching ratio: 

In contrast, if A^ were a Dirac neutrino, the measured process would correspond solely 
to the Lepton Number conserving decay vr"*" -^ e^e^fi~i'e, and the muon spectrum would 
be given by eq. (2.8), with F^v given by eq. (2.2) with C = 1, appropriate for a Dirac A^. 



The corresponding spectral and total branching ratios would then be given by: 

x^S2_M2|lii;^(m^ + m2) - ^mN{2El + ml) 



and 
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respectively. 

3 Numerical studies of branching ratios and spectra 

Here we want to analyze numerically the spectra and branching ratios for the processes 
TT^ — )• e~^e"*"^"i^ that we calculated in the previous section, to see whether they can be 
accessed by experiment within the specified range for misf and the possible values of the 
mixing elements Bin. As shown in the previous section, we express spectra and rates as 
branching ratios, relative to the charged pion decay width. 

In order to study the behavior of these quantities as a function of the neutrino mass, 
771 TV, it is convenient to factor out the dependence on the mixing elements. Here we do this 
by defining the reduced differential branching ratios as: 

dW{Tr+ -^ e+e+n-u) = .. .^ — — ^d5r(7r+ -^ e+e+^-i/), (3.1) 

{\BeNr + \Bfj.Nr) 

and similarly for the integrated branching ratios as well. We recall that /C (~ 1) is the 
correction factor defined in eq. (2.2), which represents possible neutral current channels in 
the massive neutrino decay width F^r. 

With this factor extracted, the order-of-magnitude size of our reduced branching ratios 
is given mainly by the mass of the intermediate neutrino itin, and not so much by the mixing 
elements Bin: the reduced branching ratios thus obtained from eqs. (2.10 -2.13) depend 
on the mixing elements at most in the combinations 

_ l^eAfP , . _ \B^n\ .„ „^ 

\BeN\ + \B^n\ \BeN\ + \Bfj.N\ 

quantities which are bounded between and 1. In Table 1 and figure 4 we present numer- 
ical values for the reduced branching ratio Br^ '[-k^ — ?• e'^e'^ ii~u), related to eq. (2.11), 
for various values of the neutrino mass rriN. The true branching ratios mediated by a 
Majorana neutrino N are obtained from the reduced ones after multiplying them by the 
factor (|i?eArp + \B^n\^)/1C, while if N is of Dirac type the required factor is 2\B^n\^ /K,. 
From the above we see that the measured branching ratio Br{'K~^ — )> e'^e'^^~v) is 
predicted to give in general similar values whether it is mediated by a Majorana or a Dirac 
neutrino. Consequently, the experimental measurement of this branching ratio cannot 
represent a method for distinguishing between the Majorana and the Dirac character of 



Table 1. Reduced branching ratio for the indicated process, induced by a Majorana neutrino of 
mass TTijv [c/. cqs. (2.11) and (3.1)]. 
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the intermediate neutrino, if this one exists in the relevant mass range. On the other hand, 
the measurement of the muon energy spectrum, i.e. eqs. (2.10) and (2.12), may allow us 
to distinguish between the two cases. 

In figures 5 we present the reduced spectra dBr^ ' /dE^ as a function of the muon 
energy in the intermediate neutrino N rest frame. In each graph, four of the five curves 
correspond to different values of the "admixture parameter" (om = 1, 0.8, 0.5 and 0.2) in 
the Majorana A^ case, and the fifth curve correspond to the Dirac N case, dBA^' /dE^. 
Each graph corresponds to a given representative value of the intermediate neutrino mass 
{tun = 0.112, 0.119, 0.126 and 0.133 GeV in figures 5(a)-(d), respectively). Note that 
miy must be between rra^ « 0.105 GeV and m^ ss 0.140 GeV for this process to occur. 
Consequently the kinetic energy of the muon in this frame runs over a very limited range, 

^' /2mN- This upper end is at most ~ 4.3 MeV if 



m. 



m 



V-- 



namely from zero up to {mi\f - 
rriN —7- m^r and vanishes as rriN 

Concerning these spectra, as already mentioned, E^ is the muon energy in the rest 
frame of the intermediate neutrino N . One may consider more practical to use the muon 
energy in the pion rest frame, which here we denote as E' (in the Appendix we treat in 
detail the transformation dF/dE'^ i-^ dT/dE'^). In figures 6(a)-(d), we show the spectra of 
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Figure 4. The reduced branching ratio Br^^''{-K^ — >■ e^e^pi v), eqs. (2.11) and (3.1), as a function of 
the neutrino mass ttijv. 
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Figure 5. The reduced difTerential branching ratio dBrjdE^ as a function of the muon energy in the 
neutrino TV rest frame, iJ^j, as defined via eqs. (2.10)-(3.1), for the decays 7r+ — >■ e^e^ ^iTv mediated by 
a Majorana neutrino N , for various neutrino masses: (a) mjv = 0.112 GeV; (b) ttiat = 0.119 GeV; (c) 
rriN = 0.126 GeV; (d) viin ~ 0.133 GeV. In each graph there are five curves, corresponding to different 
values of the admixture parameter aM [eq. (3.2)]: om — 1.0 is the soHd (M) curve; 0.8 (dashed); 0.5 (dot- 
dashed); 0.2 (dotted). The case mediated by a Dirac neutrino is also presented as the soUd hue labelled 
(D), with the distribution normalized so that the area under the curve is the same as in the Majorana cases. 
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Figure 6. Same spectral distributions as in figures 5, but now for the muon energy S^ measured in the 
rest frame of the pion. 



the same cases, respectively, but as functions of the muon energy in the pion rest frame, 
dB^/dE'^. 
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The om = 1 case is denoted in these figures by the sohd fine with the label (M), be- 
cause for this admixture the Lepton Number conserving process tt"*" -^ e^e~^fj,~iye vanishes, 
and so the only contribution comes from the Lepton Number violating process, mediated 
necessarily by a Major ana A^. Complementary, the purely Dirac neutrino N case is shown 
as the solid line with the label (D). 

These figures suggest that measurements of the differential decay width dT/dE'^ can 
more or less distinguish the Majorana vs. Dirac character of the intermediate neutrino 
N, the distinction being clearer when the "admixture parameter" aM approaches unity 
(i.e. l-BeAfP ^ l-B/iATp, so that the Lepton Number conserving component Tr"*" -^ e+e"'~/u~z^e 
is relatively suppressed in the Majorana case). Furthermore, comparison of the results 
in figures 5 with the corresponding results in figures 6 shows that the difference can be 
discerned apparently more clearly in the N rest frame than in the pion rest frame, provided 
the muon energy can be measured with enough precision: in the A^ rest frame the curves 
for the purely Dirac cases are increasing with energy Ef^ all the way to the upper endpoint, 
where they suddenly drop to zero (the sharp drop is smoothed out by the finite electron 
mass), while the curves for the Majorana cases reach earlier a maximum and then drop 
gradually towards the upper endpoint. In contrast, in the pion rest frame the muon spectra 
do not show such a clear distinction between the Dirac and Majorana cases. 

Now, concerning the experimental challenges, the determination of the muon energy 
in the lab frame needs to have an uncertainty below 1 MeV to achieve the required level of 
discrimination. This requirement is realistic with current detector technology (for example 
the momentum resolution for a muon in the inner detector of ATLAS is near lO"'^ [40], 
which means that a 1 GeV muon can be measured with a precision of a few times 0.1 MeV). 
Now, in the pion rest frame or in the A^ rest frame, the muon energy, E' and E'^ respectively, 
is obtained after further kinematic analysis. In general the former may be obtained with 
more precision, because it can be inferred just from the measured muon momentum in the 
lab frame, provided an accurate energy of the pion beam is known; instead, for a precise 
determination of the muon energy in the neutrino A^ frame, a precise measurement of the 
primary positron momentum is needed as well. In any case the required precision for the 
positron momentum is also realistic for today detectors. Probably the most challenging 
issue here is to have a pion beam with an energy spread, 6Et^/Et^ (in the lab frame) which 
is small enough. For highly relativistic pions, the corresponding uncertainty of the muon 
energy in the pion rest frame is 5E' /E' Ki SEt^/Et^. This uncertainty needs to be < 10"^, 
so that the energies E' can be determined with an uncertainty below 1 MeV, as required. 

In future high intensity beam facilities, such as Project X at Fermilab, charged pion 
beams with lab energies Et^ '--'2 — 15 GeV and a luminosity '--' 10^^ cm~^ s~^ will be 
produced [41]. Taking a beam of 1 cm^, we can expect '^ 10^^ pions per day. With such a 
large sample of pion decays, the search for the primary process tt~^ — > e^A^ [cf. eq. (2.6)] 
can already improve the current bound on the mixing element |i?e/vfp ^ 10^^ [39] by many 
orders of magnitude. 

According to Table 1, for |i?eArp ~ 10~® we get branching ratios Br*^ •'(vr"^ -^ e^e^ ^~v) 
10~^^ when the neutrino mass mjv is in the range between 0.120 and 0.135 GeV, and thus, 
in principle, 10^^ such events could be detected per year. However, as explained in the 
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previous section, the probability of the neutrino N to decay inside the detector, with the 
detector of length L ~ 10^ m, is 

Pn ~ L/{-f^CTN) ~ l^eTvl' ~ (lO"" " IQ-^^) l^eiV |' , (3.3) 

where 77r ~ 10^ — 10^ for E^^ ^ 2 — 15 GeV. This acceptance factor, for |-BeArP ~ 10~^, 
is thus expected to be ~ 10"^^ — 10"^^, meaning that, instead of the ~ 10^^ events just 
mentioned, about 10^ — 10^ such events can be detected per year. Now, if |i?eArp is smaller 
by a factor '^ 10, the number of such detected events would be lower by a factor ~ 10^, 
because Br oc |i?eAfp and the acceptance P/v oc |i?eAfp5 cf. eqs. (2.10)-(2.11) and (3.3). 
Keeping in mind the current upper bounds for jSeTvP ^ 10~^, we see that the search for 
the events of the type vr^ — )• e~^e~^ fi^i' is promising, due to the very large expected number 
of produced charged pions. If |-Be^P ^ 10"^*^, the number of such events would be > 10^ 
per year, which may allow us to distinguish between the Majorana and Dirac character of 
neutrino N, cf. figures 5 and 6. 

4 Summary and Conclusions 

We have studied the possibility to discover a sterile neutrino A^ in the mass range between 
771^ and rriTr, and detect its Majorana or Dirac nature, using the charged pion decays 
tt'^ -^ e'^e^ ^"v^ and vr"*" -^ e^e^fi~i'e- The neutrino in question is in the intermediate 
state, and in order to have sizable experimental signals, it must be on its mass shell (hence 
the specified mass range). Even so, the rates are extremely small, so they can only be 
detected in extremely high intensity pion beam experiments. The first process violates 
Lepton Number by two units, so it can only be produced if the intermediate neutrino N is 
Majorana. In contrast, the second process only violates Lepton Flavor, while conserving 
Lepton Number, so it can be produced indistinctly by a Majorana or Dirac neutrino N. 
However, given that the final neutrino fiavor (i>^ or i/g) is not experimentally observed, 
both processes can contribute to the measured signal, which we refer to as vr"*" — > e~^ e~^ fj,~ u . 
Moreover, the rates of both processes could be comparable, so in general the measurement 
of the branching ratio Br{'K~^ — >• e"'"e"'"/i~z^) allows the discovery of the neutrino N , but 
does not distinguish its Majorana or Dirac nature. 

While the rates vanish if the neutrino mass mj^j is at one of the ends of the range (w.^, 
m-,^), for 771 TV near the middle of this range the branching ratio can be up to Bv^tt'^ — > 
e'^ e~^ fi^ v) ~ |i?£Arp X 10~^. For the current upper bound on the mixing, namely I-B^at^ < 
10^^, this means that the branching ratio can be up to 10~^^. 

Since a neutrino N of the required mass must be sterile, its lifetime is rather long, so 
the decay separates in time (space) as tt"*" — > e"'"A^ followed hy N ^ e~^ fi~ u . Consequently, 
the first signal for the discovery of such a neutrino in these kinds of experiments comes 
simply from the detection of the prompt positron with an energy well below the standard 
mode vr"^ — > C^Ue in the pion rest frame. 

While the decay of this standard mode is chirally suppressed due to the small positron 
mass, the non-standard primary process vr"*" — >■ e"'"A'^ is only suppressed by the mixing 
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element B^n (actually looking for low energy bumps in the positron spectrum is how the 
current bound on this mixing element was obtained). Consequently, the absence of a low 
energy bump in the primary positron energy spectrum in these high intensity pion beam 
experiments will allow to improve by many orders of magnitude the existing upper bounds 
on the mixing element |i?eAf| (currently |i?eArp < 10^'^, for mj\j in our range of interest). 

The Majorana or Dirac nature of the neutrino N cannot be determined from the 
branching ratios or from the pure detection of the prompt positron. We have studied the 
possibility to determine such feature from the energy distribution of the muon (a muon 
with charge opposite to that of the decaying pion). We find that the muon spectrum 
in the rest frame of the intermediate neutrino N is quite different for the two processes 
7r+ — )■ e'^e^fi^L'f^ and vr'^ -^ e~^ e'^ fi~ Ue'- while the first reaches a maximum and then 
decreases gradually to zero at the upper endpoint, the second grows monotonically all the 
way up to the endpoint, where it sharply drops to zero. If iV is a Dirac neutrino, the 
second feature will clearly show, as the first process will be forbidden (a Dirac neutrino 
cannot induce Lepton Number violation). On the other hand, if A^ is a Majorana neutrino, 
both processes will occur and the distinction could be less clear, depending on the mixing 
parameters. If {B^nI ^ l-^^A^I (i-e- OiM — ^ 1)> the Lepton Number violating process 
TT^ -^ e~^ e~^ fi~ Ufj^ dominates, and the spectrum will show more clearly its shape, signaling 
the presence of a Majorana neutrino. On the other hand, if l-BeAri ^ {B/j^n], it is the Lepton 
Number conserving process that dominates, so that even if A^ is Majorana, the spectrum 
will show the same shape as if N were a Dirac neutrino. In either case, a good energy 
resolution will be important to tackle the separation between the Dirac and Majorana 
nature of the neutrino N using the muon spectrum. 
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A Appendix: general formulas for M+ — ;■ ifN — )• i'^^J^"// in the Af+ and 
A^ rest frames 

In this Appendix we present general formulas for the the Lepton Number violating (LNV) 
and the Lepton Number conserving (LNC) decay of a charged pseudoscalar meson M+, 
M+(pAf) — >■ ■^i'(pi)^2'(^2)-^~(j'^)^(Pj/); cf- figures 1 and 2, respectively. Both decays are 
assumed to take place via the exchange of an on-shell neutrino N which is assumed to be, 
in general, Majorana. The masses of the particles are assumed in general to be nonzero, 
exccept for the mass of the standard neutrino (rrijy = 0) . For the direct channel of the LNV 
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decay M~^ — >■ ^^^^^ ve (cf. figure 1) the transition amplitude is 

-^ = 2 T-—- r K(P^)7^(1 - l^)vv{Vu)\ ^2(^2)7^24/(1 - l5)viAPi)\ 

Pn ^^n ' iTni^i jv 

(A.l) 
We use the notations: T^ is the total decay width of the neutrino A^; Gp is the Fermi 
coupling constant {Gp ~ 1.166 • 10~^ GeV~ ); -B^ at is the mixing element between the 
neutrino of flavor state vi- with the (mass eigenstate) neutrino N , cf. eq. (2.4); X^ is the 
phase factor of the Major ana neutrino N {\\n\ = l)j fhi is the decay constant of the meson 
M"*"; and VqQ is the Cabibbo-Kobayashi-Maskawa (CKM) matrix element corresponding 
to M"*". The subscript c in the spinor for £^ in eq. (A.l) means the charge-conjugation 
[v^ = u m. the helicity basis). Furthermore, we will denote by rrij the mass of t^ lepton 
(J = l,2); 

Since the process is dominated by the intermediate neutrino A^ on mass shell, its prop- 
agator in the transition probability |7Wp can be given by the narrow width approximation: 

^ ^ 5{pl-ml). (A.2) 



For the calculation of the decay rate, T{M~^ — > £^£^£~i>£) = {2mM)~^ f dps4\A4\'^, it is 
convenient to decompose the 4-body phase space integral as the 2-body phase space of pi 
and Pn (where the intermediate momentum p^ has a variable mass hn), and the 3-body 
phase space of p2, Pi and p^ following from the A^ decay: 

d/iAT^ f Spi d?PN ,„ .4.4. ^ .. „s 

{2Tiy5*{pM -Pi-PN) (A.3) 



27r J {2-Kf2Ei {2T^f2EN 

d?p2 (fpe d?py , .44, . 

i2'nnE,{2.Y2E,{2.Y2E} ^''^ ^ (PiV -P2 -p, -p.), 

where the range of the /uat integration is [i^ £ ("i^ + "i2>"ij\/ — mi). This is precisely 
the range of possible neutrino masses tti^t where A^ can be on mass shell in the process. 
We can now integrate the probability over the phase space. Since the integral expression 
is Lorentz invariant, we can choose any reference frame. It turns out that, in the narrow 
width approximation (A.2), the direct channel (figure 1) of the LNV process and the crossed 
process (pi o pi and m\ -H- 7712) contribute as a sum to T, i.e., the interference term is 
zero. 

We work, for convenience, in the center-of-momentum (rest) frame of A^, CM(A^), 
where pjq = (7nAr,0).^ We choose the z direction along pi, denote the angle between pi 
and Pn as 9^, and the angle between pi and pi as 61^ see figure 7 (we refer to the diagram 
of the type of figure 1). For the LNV decay M"*" — t- ^^^^£~^'£ (cf. figure 1 and the crossed 



■^ Nonetheless, the (differential) decay widths dT will be written for the M"*" rest frame, i.e., without the 
time dilation factor. 
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Figure 7. 



diagram) we obtain (in CM(A^)) 



dEid cos 9 i 



^{LNV) 



[rriNimN - 2En) + 



nil — m 



212 



4:7nj\[ [mj\f{mN — 2E() + 



nil 



X <^ cosdi{m% - ml) [(m|^ - m%)'^ - 2m\{m\.i + ni%) + mf\ {EJ - nij) 
+ [m%{nilj - ni%) + m\{ni\j + 2m%) - nv{\ Ei>\pi>\ \ + (mi o 771-2) , (A. 4) 



where Z^ ' is defined as 

and A^' ^ is the square root of the function 

\{x,y,z)=x +y +z —2xy — 2yz — 2zx 



m'ij,ni%,ml) , (A. 5) 



(A.6) 



We note that the factor (1 — 61-^^12/2) in Z^^^^' accounts for the fact that, when if 7^ l^, 
there are two types of decays (with intermediate on-sheh A^) leading to the result M+ — )• 
ifip-iyf. (a) M+ -^ i+N -^ i+£+i-i?i, and (b) M+ -^ e+N -^ £+£+£-9^. 
Integration over the angle Oi gives 



dEp 



(M+ ^ et^rue) = ^^^^^^^ [mUm% + ml) - {m% - ml)'] 



{m%-2mNEe + mj-miy 
xi^f\ Hji —m„ 



(m^ - 2mMEi + m|) 



+ (mi f^ m2) , 



(A.7) 



where the first ("d") channel term (for: M^ -^ ifN — ?> £^£2^ Vf) contributes when 
me,<Ei,< (£^£)i^L, and the second ("c") channel term (for: M+ -^ £+N -^ £+£+£' i?e) 
contributes when m^ < Ei < (£'£)max, where 



{Ei 



)(d) 
/max 



2m, AT 



/' 2 I 2 2N 

(m^ + m£ - mgj 



(^, 



/max 



2m. JV 



(m^ + 771^ — 7ni) . (Ai 
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When the masses mi and m2 are neghgible, we obtain 



dE id COS 9 

dYiLNV) 



dEi 



Z^^^^^ -mN{rn\i-m%)\pi\ [mj + mN{mM - 2£'^)] [cos Oe\pe\ + Ef] , 

(A.9) 
Z^^^^^mN{m\i-m\)Ei^ Ej - mj {m% - 2mNEi + m|) , (A.IO) 

r(-^)(M+ -. Wr..)L^_^, = z(-^)^ (^ - l) /KVm^) , (A.ll) 

where f{x) is the 3-body decay function, eq. (2.1). 

For the LNC decay M+ -^ t[l2^~Uj^ (cf. figure 2 and the crossed diagram) we obtain 



mi=m2=0 



mi=m2=0 



iE(.d COS Oi 24mN [mj + mN{mN - 2Ee)f I 



m%) 



dEed cos 0f 24mN [mj + mN{mN - 2Ei) 

cos9i>{m\ - m\!)\pi\^J {{mM + miY - m%){(mM - miY 

X [ (37TT,£ + m]\[{m,]\[ — 4£'£)) (jn^ + m,]\[{m,N — 2Ei)) + m,2 (Sm,^ — m,]^{m.]^ + 2Ei)) ] 

+ [ (mf - m^(mi, - m^) - ml{ml, + 2m^)) {SEJml - 2mjmr,{2ml + mj + m%) 

+2E^mN{m2 + 5m^ + 5mj^) + E'^(3?77.2^ + 3m2m^ + (3Tn£ + m,j>^){m^ + 3mjv)))] 



+(mi om2,4 ^^2) , (A.12) 

where Z^ ' is defined as 

ziLNC) ^ GUBl^B,r^V:Q\^ (1 - i^,,,,) (^^AV2K,,m^,m?) . (A.13) 
Integration over 61 gives 



dEe '"^ .^i^^^-/ 96m^ [^2^m^(_2^^ + m^)]^ 

X < 8w (E'l — m,j)mN [m\ — m,j + (2£'^ — ?7T,Ar)7njv] 

X [— mf + m\^m,j^ — m,% + m\{mi\^ + 2m^)] [8£^|m^ — 2m|m7v(2m| + m| + m^) 
—2E(^m,N (^-2 + 5(m£ + rn^v)) + Ei{3m,£ + lOm^m^v + 3?7ijv + 3m2(r?T,£ + m,^))] > 
+(4 ^^2, "11 07712) , (A.14) 

where the exphcitly written ("d") channel term (for: Af+ — > ^^A^ o- if £~ £2 1^2) has 
nonzero values in the interval mc < Ei < (Sf)max, and the second ("c") channel term (for: 
M+ -> ^^A^ -^ ifl'ifui) in the interval me < Ei < (£'^)max, where (£'^)max (2; = d,c) are 
given in eq. (A. 8). 



^the cases j = 1,2 were added 
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When the masses mi and ?7i2 are neghgible, we obtain 



driLNC) 



dEfd cos 6f 



mi=m2=0 



ZiLNC) ^ (^^ ^ ^^)^ ^m^iml, - m%M\ 

X < — 3m^ (cos 6*^1^^! + 2mN — Ei) + 

mN [cos ee\pe\{4:Ee - iun) + {-^Pi\^ + Sm^v^^)] } , (A.15) 



driLNC) 



dEf 



mi=m2=0 



(^Z(^NC) ^ (^^ ^ ^^)^ \mM{mlj - ml)^EJ-MJ 



{m% + mj)Ei - -niNiEJ + mj) 



(A.16) 



r(--)(M- ^ iUtr.)U.^^^, = \ (z(--) + (ii o ..)) ^ (^ - 1) /KV-^), 



(A.17) 

the last expression being almost identical with eq. (A. 11) for the LNV decays, except that 
now Z(^^^) of eq. (A.5) is replaced by Z^^^^) of eq. (A.13). 

The general case of Majorana (M) neutrino N includes LNV and LNC decays, and in 
such a case dr(*^) = dV^^^^~^ + dT^^^^~^ (if £i / £ and ^27^^)- 

We note that the case of negligible masses mi and ?7i2 and with (.1 = ^2 is applicable 
to the decays considered in the main text of this work, namely £^ = £^ = e"^; £~ = /x~; 
and M^ = vr"*". Specifically, the expressions (A.lO)-(A.ll) and (A.16)-(A.17) acquire in 
this case the following explicit form: 



rfvr-' 



e^e^/i v^) 



dE„ 



(A.18) 



G\\B^fl iKdP m% {ml - m%f 



167r^ ?7i^ Ftv 
_G^|Se^|V.'|14dpm^(m2-m^)2 m% 



Ef, {m% + ml- 2mNE^)JEl - ml, 



167r4 ml V 



N 



96 



f [ml/ml) , (A.19) 



li IR ,,R* \2f2 IT/ .12^2 



m'j^+rn^ 



dE„ 



(A.20) 



2 \2 



G^|i?eivi?;^|^/^ iKdl' m^(m^ - mj,) 



32tt^ ml r 



AT 



2 
— ? 
3 



(m^ + ml)E^ - -mN{2E'l + m^' 



El - ml, 



p. +^ + _ + ^ G%\B,NB*,,\''fl \Vud? m%{ml - m%f m% , ^ 2^ .. .^^ 
r(vr ^e ;. e ..) = isvr^ rn3 r^ ^ / K^^) ' (^-21) 

where /(x) is given in eq. (2.1) and Tat is given in eq. (2.2). 

The differential distributions dT/dE^ for the decay M+ -^ if £2 '^~ i'i{i'j) , eqs. (A. 7) and 
(A. 14), refer to the quantities Ei in the center of momentum (rest) frame of the intermediate 



17 



neutrino N: T, = CM(A^), see also figure 7, while the dT refers to the differential of the 
decay width in the center of momentum frame (rest) of the decaying particle A/"*": S = 
CM(M+). We define the reference axis z in both frames as the direction of pi of if, for 
the decays of the type of figures 1, 2. 

The initial meson M"*" is without spin. Therefore, for the direct channel (cf. figures 
1 and 2) the distribution dr/{dE^d9 ) in the frame S = CM(M"'") is independent of the 
angle 6 , where 9 is the angle between a chosen reference axis z and the direction of if in 
CM(M+) 

E=CM(N) (a) ^ =CM(M^ (b) 

^ .N> Pi(-Pm) , ,,.^ ? ,,,, . Pn .. ^'""\ ,.. ? 
x-t:^ :r- -f/iJ (N)- -^>^ ^ — - (/, j 



e, ; 



(M) ^.0'^ 



p/ ~~-fn P' \f/"j 



Figure 8. 

dr 1 dr 



— f A 22) 

dE'^dn{6',(l)') AirdEy ^ ' ^ 

Therefore, we can choose this angle to be ^ = 0, i.e., Pi = \pi\z, cf. figure 8(b). This axis 
z is the same also in CM(iV), figure 8(a), i.e., z = pi = Pi- 

In E = CM(M"'"), the momentum p^ in the considered direct channel is 

E'n = i^i^i'nnld + ml, - ml) , p'^ = -|p^| = -^-— A^/^(mi^,m^,m?) .(A.23) 

The components of p^ in the two frames S = CM(A^) and S = CM(M+) are then related 
via a boost determined by p^. The angle between pi and pi (= z) is ^^ in S = CM(A^), 
and 0^ in S = CM(M"''), i.e., p| = \pe,\ cos Oi and p^ = \p^\ cos 6^. Then, E^ and cos6'£ can 
be expressed as functions of E^, cosO^, Ej^ and p^, where the last two quantities are fixed, 
cf. eqs. (A.23). In this way, we obtain the needed one-to-one correspondence between the 
energy and the (azimuthal) angle of l~ lepton in the two systems CM(A^) and CM(M"'") 

Ei = Ee (^E'e,cose'^^ = -^ (e'^e'^ - Pn yJiK^ - mj) cos O'A , (A.24) 



cos 6i = cos Oi (E^, cos 



-PnE^ + E^\pi\cos6^ 



e'nK - Pn \J{K^ - mj) cos e'A 



2 

2 2 

— mj^rrii 



-1/2 



(A.25) 



The distributions dT/dEi and dT /dE^ are then related via the following relation: 



dcosOi—-— —\J\ , (A.26) 



^£ 



dE'p Jc' . dEid cos i 
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where c^^^ = [cos{9^{Eg))]jnm, Cmax = [cos(0^(ii^^))]max5 and J is the corresponding Jacobian 

_ djEj, cos 9i) _ dEj d cos 9^ dEj d cos Oj f A 271 

d{E'f, , cos 9'^) dE'^ d cos 0^ d cos 6*^ 9£'^ 

which can be calculated from the obtained relations (A.24)-(A.25). The result for J is 
J = ruN Upel [e'^e'^ - mjinj^j - 2E[E'j^\pf\'^p^cos9[ + \p(\^{E'^ - m^)cos^6'H 

-rnlml, + (e'^e'j^ - \p^\p^cos9'^ \ (A.28) 



We recall that E^ and p^ are fixed, eqs. (A. 23). 

The actual differential decay width dT/dEf^ as measured in the T, = CM(M^) frame 
is then obtained by integrating the expression (A. 26) over all the angles 9^, where 9^ is the 
angle (in S ) between if and £~ for the decays of the type of figures 1, 2, see figure 8(b). The 
expression to be integrated over 9^ is the double differential decay width dT /{dE^d cos 9 n), 
where En and cos^^ are now considered functions of E'^ and cos9g, eqs. (A.24)-(A.25). The 
double differential decay widths entering the integral in eq. (A. 26) are given in eqs. (A. 4) 
and (A. 9) for the LNV decays, and in eqs. (A. 12) and (A. 15) for the LNC decays. Since 
the relevant parameters in the integration in eq. (A. 26) are now E^ and cos^^, we also need 
the relations (A.24)-(A.25). 

The integration limits in the integral (A. 26) can be obtained in the following way. The 
relations inverse to (A.24)-(A.25) are 

e'^ = (E'j^Ee + p'^\pe\ cos 9e) = E'^{Ee, cos 9'^) , (A.29) 



niN 



'3 



COS 6*^ = [p^Ee + E]sf\pe\ cost 



(e'^Ei + Pn\pi\ COS 9 A - mj^mj 



-1/2 



. (A.30) 



Prom eqs. (A.29)-(A.30) we obtain for dcos9^/dEi at fixed E^ a positive expression, 
"i^N/lPel/i—PN)- This implies that Ei grows when cos 9^ grows, at fixed Eg. This, in 
conjunction with eq. (A. 24), means that, at fixed E^, the quantity rriNEi varies between 
{Ej^Eg+p^\pg\) and {Ej^Ef.—p^\pg\). Simultaneously, we know that Ei must vary between 
m£ and (ii^f)max of eq. (A. 8) (in the direct channel). Therefore, Eg, at given fixed E^, varies 
within and covers the interval which is the overlap of the two aforementioned intervals, i.e., 
Ei{E[)rain < Ef < Ei{E'g)mi,^, where 

mNEi{E'g)rnin = Max(mNmi,E',^E'g + p'^\pA , (A. 31) 

mNEi{E'()ra-s.^ = Min i -{mj^ + mj - ml),E'^E'g - PnIPiI j ■ (A.32) 

To obtain the corresponding minimal and maximal values of cos9g, at given fixed E^, we 
use an expression of cos 9^ in terms of E^ and Ei, which is obtained from eqs. (A.29)-(A.30) 

cos 6*^ = , =- — TT-imNEe - E^E^) = cos9g{Eg, En) . (A. 33) 

Ie7 - m? (-Pn) 
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Accounting for the aforementioned fact that, at fixed E^, the quantity cos0£ grows when 
El grows, we finally obtain the following minimal and maximal values of cos^^, at given 
fixed E^, i.e., the lower and the upper bounds of integration in eq. (A. 26): 

Cmin = C0S6'^ ( E'^, Ei{E'i),am) = , =7 7n7 ( mAr£'£(£^)mm " E'jsfEA ,(A.34) 

^ ^ ^E'/-mj^-PN)^ ' 

cinax = COS 6*^ [e!i,Ei{E!()^^A = ^ ( mAr£^K-S£)max - E'j^E'A (A.35) 

^ ^ ^42 - mj ^-Pn) ^ ' 

where the energies Ei{E^^-^^ and Ei{E^)^ax, for a given value of E^, are given in eqs. (A.31)- 
(A.32). This specifies the integration limits in the integral (A. 26). 

We note that in the frame S = CM(M"'") the energy of the i~ lepton E^ varies between 
m£ and the energy when Eg is maximal [eq. (A. 8)] and cosO^ = — 1 



nil ^ Ep < 



2m% 



E'j^{m% + mj-ml) + {-p'^)X^/^{m%,mj,ml)\ . (A.36) 



Using all these formulas, we can perform explicit numerical evaluations for the decay 
M+ -^ £f£2^~i'f for the case of the pion decay into two positrons: M+ = 7r+, £~ = ;U~, 
£^ = £2= c"*" (^-i = "22 = 0). Here, Uf = v^^Vf. when the intermediate on-shell neutrino 
N is Majorana and Dirac, respectively. The mass rriM, for on-shellness, is required to 
be between m^ = 0.10566 GeV and 771,^ = 0.13957 GeV. Details of the results of such 
calculations are given in section 3. We checked numerically that the integration of the 
differential decay width over the energy of the ^~ lepton gave us the same result in the 
CM(7r"'") as in the CM(A^). This represents a strong cross-check that our formulas for 
calculation of the differential decay width dV^^' /dE^ in the meson rest frame CM{M~^) 
[=CM(7r"'")], i.e., eq. (A. 26) and the subsequent formulas, are correct. 
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